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Objective Decreased heart rate variability (HRV) has been reported in generalized anxiety disorder (GAD), but the results are mixed. 
Little is known about the impact of comorbid major depression (MD) on HRV in GAD patients. Both issues necessitate further investigation. 

Methods Twenty unmedicated, physically healthy GAD patients, 20 GAD patients with a secondary diagnosis of MD, 40 MD patients 
and 60 matched controls were recruited. We used the Hamilton Anxiety Rating Scale and the Hamilton Depression Rating Scale to as- 
sess anxiety and depression severity, respectively. Cardiac autonomic function was evaluated by measuring HRV parameters. Frequency- 
domain indices of HRV were obtained. 

Results Three patient groups had more anxiety and depression symptoms than control subjects, but heart rates (HRs) were signifi- 
cantly elevated only in GAD patients with comorbid depression. Relative to controls, GAD patients had reduced HRV while GAD pa- 
tients with comorbid depression displayed the greatest reductions in HRV among three patients groups. Correlation analyses revealed 
anxiety/depression severity significantly associated with HRs, variance, LF-HRV and HF-HRV However, separately analyzing among 
individual groups and adjusting for HRV-associated covariables rendered the correlations non-significant. 

Conclusion Our results suggest that reduction in HRV is a psychophysiological marker of GAD and individuals with comorbid GAD 
and MD may be distinguished based on psychophysiological correlates (for example, HF-HRV) from non-comorbid GAD patients. 
Taken into account that comorbid depression may confer increased risks for cardiovascular events in GAD patients, this subgroup of 
GAD patients may benefit better from cardiovascular risk reduction strategies. Psychiatry Investig 201 3; 1 0:326-335 
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INTRODUCTION 

Generalized anxiety disorder (GAD) is a chronic anxiety 
disorder and characterized by fluctuating levels of persistent, 
uncontrollable worry associated with fatigue, insomnia, mus- 
cle tension, poor concentration, and irritability. Anxiety dis- 
orders have been associated with increased risks of cardiovas- 
cular morbidity. 1 One potential explanation proposed for this 
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association is a dysregulation of the autonomic control of the 
heart because autonomic nervous system (ANS) activity is as- 
sociated with cardiovascular disease (CVD) 2,3 and anxiety 
disorders. 1 Many previous studies has examined the impact of 
phobic anxiety on patients with established coronary artery 
disease (CAD), 1,4,5 few examine GAD-specific contribution. 
Until recently, there were two published GAD-specific studies 
quantified GAD as having a distinctly negative impact on 
outcome of patients with stable CAD. 6,7 It raises the question 
of whether the relationship between GAD and the reported 
increase of cardiac risks is mediated by cardiac autonomic 
dysregulation. Because CVD itself might influence cardiac au- 
tonomic function, we focused on GAD patients without CVD, 
in order to avoid overestimation of the association between 
GAD and cardiac autonomic function. 

The ANS links the central nervous system and the cardio- 
vascular system. Efferent links in the neural control of heart 
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rhythm consist of sympathetic and parasympathetic fibers in- 
nervating the sinus node. Because sympathetic and parasym- 
pathetic firing alters spontaneous sinus node depolarization, 
cardiac rate and rhythm convey information about autonomic 
influences to the heart. 8 Previous research focused primarily 
on the sympathetic branch of the ANS and reported higher 
levels of sympathetic activity in GAD patients. 9 " 13 However, 
few psychophysiologic studies of GAD use a technique en- 
abling quantitative assessment of influences of the parasym- 
pathetic nervous system (PNS), the other branch the ANS. 
The PNS controls heart rate (HR) by brain stem regions that 
project to the heart via the vagus nerve. These vagal fibers 
functionally slow HR and actively inhibit the sympathetic in- 
fluences to the heart. 14 Vagal nerve activity thus leads to regu- 
lar oscillations in HR Cardiac vagal control (CVC), essential- 
ly, reflects the degree to which there is tonic vagal influence 
on the heart. 15 Low CVC, which causes reduced autonomic 
flexibility, has been associated with a number of psychopatho- 
logical states, including depression, 15 anxiety, 16 and the cogni- 
tive processes that have been associated with depression and 
anxiety, namely worrying. 17 Given the central role of worry in 
the DSM-IV diagnostic criteria for GAD, CVC has recently 
become a promising candidate for investigation. Our study 
used frequency-domain analysis of heart rate variability (HRV), 
a method different from previous studies, 9 13 providing indi- 
ces that reflect sympathetic and parasympathetic cardiovas- 
cular influences. 

Frequency-domain analysis of HRV, with its standard pro- 
cedure and interpretation first reported in 1996, is a sophisti- 
cated and noninvasive tool for detection of ANS regulation of 
the heart. 18 An important advantage of frequency-domain 
analysis of HRV is that it utilizes spontaneous fluctuations in 
HR to estimate tonic ANS functions. However, it should be 
noted that a well controlled condition is required for sponta- 
neous ANS functional recordings. In the present study, we 
used spectral analysis of HRV while subjects were in supine 
rest in a quiet and relaxing atmosphere for 5-min short re- 
cording times. We did this because short-term resting-state 
measures of HRV better reflect intrinsic HRV. 19 

Studies have seldom examined the impact of GAD on HRV. 
As far as we know, the results of those studies concerning rest- 
ing HRV in GAD are mixed. For example, some studies found 
GAD patients exhibited an increased HR and decreased car- 
diac vagal tone during rest as compared to controls. 17,20,21 How- 
ever, Licht et al. 22 reported that HRV reductions in GAD pa- 
tients are driven by medication effects alone rather than the 
diagnosis of GAD. The conflicting results might be due to 
heterogeneity in relatively small samples, confounds from 
medication, physical health, habitual physical activity, smok- 
ing, psychiatric comorbidities, and reporting of different HRV 



measures. To address these concerns as mentioned above, one 
should analyze with better methodology. GAD is often asso- 
ciated with an increased risk of developing a secondary diag- 
nosis of major depression (MD), as suggested by epidemio- 
logical studies. 23,24 For example, around two thirds of patients 
with a lifetime diagnosis of GAD retrospectively report MD, 25 
suggesting that GAD has a tendency to precede depression 
and finally develops into depression. However, decreased 
HRV can be seen in MD 26 and the measures of HRV have 
been correlated with depressive severity, including both glob- 
al 26 and symptom-specific indices. 26,27 Thus, factors concomi- 
tant with the MD diagnosis itself might influence the observed 
relationship between GAD and HRV. One might wonder 
whether the comorbid diagnosis of MD in GAD patients is as- 
sociated with greater reduction in HRV as compared to GAD 
patients without comorbidity. Taken together, the following 
hypotheses were tested: 1) Physically healthy, unmedicated 
GAD patients without comorbidity will have decreased rest- 
ing HRV as compared with age- and sex-matched controls; 
and 2) GAD with comorbid MD will display greater reduc- 
tions in resting HRV as compared to non-comorbid GAD/ 
orMD. 

METHODS 
Subjects 

This study was approved by the Institutional Review Board 
for the Protection of Human Subjects at the Tri-Service Gen- 
eral Hospital, a medical teaching hospital of the National De- 
fense Medical Center in Taipei, Taiwan. We obtained written 
informed consent from all participants and fully explained 
the procedures of the study. Initial study entry criteria: age 20- 
65. After detailed questionnaire screening, clinical examina- 
tion and chart review, we excluded subjects with pregnancy, 
smoking, diabetes, cancer, neuropathy, any cardiovascular 
disease that affects HRV or engaging in regular physical train- 
ing exceeding 10 hours a week. Subjects who used any medi- 
cation (e.g., antipsychotics, anticholinergics,antidepressants, 
oral contraceptives, anticonvulsants, anxiolytics, cerebral met- 
abolic activators, or cerebral vasodilators) that have been re- 
ported to affect the ANS functioning for at least two weeks be- 
fore evaluation were also excluded. 

Based on the same methodology in our previous studies, 26,28 
each patient was evaluated using the Chinese Version of the 
Modified Schedule of Affective Disorder and Schizophrenia- 
Lifetime (SADSL) 29 to reach DSM-IV criteria for a primary 
diagnosis of GAD. Our previous study has reported diagnos- 
tic data with satisfactory interrater reliability 30 Here, we fur- 
ther excluded individuals with a history of substance depen- 
dence, organic brain disease, any comorbid anxiety disorders 
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or major psychiatric disorders other than MD. GAD patients 
without comorbidity were recruited in GAD group (n=20). 
GAD patients who met criteria for an additional current di- 
agnosis of major depression (MD) were recruited in GAD 
with MD group (n=20). The severity of MD was assessed 
with the 17-item version of the Hamilton Depression Rating 
Scale (HAM-D). Only subjects with a minimum score of 18 
on the HAM-D entered the GAD with MD group. We also 
recruited the "depression-only" group from outpatient set- 
ting. This group comprised MD patients without comorbidity 
(n=40), who are parts of participants described in one of our 
previous studies. 26 Finally, control group (n=60) was recruited 
from the community. We used the modified Chinese Version 
of SADSL to exclude individuals with psychiatric conditions. 
Control subjects had no lifetime history of a mental disorder. 
Three patient groups were matched to control group on age, 
gender, and education. 

Assessment of anxiety severity 

All participants were assessed using self-report measures of 
anxiety, i.e., Beck Anxiety Inventory (BAI). 31 The BAI contains 
21 items that measure anxiety-related symptoms and was 
shown to have good test-retest reliability and fair concurrent 
validity Subjects were asked to rate the severity of their anxi- 
ety on 4-point Likert scales (0-3). They were also assessed by 
an attending psychiatrist (HAC) using clinician-rated scales, 
i.e., the Hamilton Anxiety Rating Scale (HAM-A). 32 Both the 
BAI and the HAM-A provide global indices of anxiety sever- 
ity. To avoid multiple testing of the same hypothesis the anal- 
ysis of the relationship between HRV parameters and global 
anxiety, severity was based on the HAM-A. The results were 
unchanged whether interviewer or self-reported measures of 
anxiety severity were used as an outcome. 

Assessment of depression severity 

All participants were interviewed by HAC using the 17- 
item HAM-D, an objective scale to assess severity of depres- 
sion. In addition, we used the Beck Depression Inventory 
(BDI), a 21-item questionnaire, to assess subjects' self-report- 
ed severity of depression. 33 The results are scored by summing 
the responses to each of the items in order to obtain a total de- 
pression score (range, 0-63). 

Measurements of HRV 

Detailed procedures have been reported previously. 34 In 
short, after sitting quietly for 20 min, a lead I electrocardio- 
gram was taken for 5 min while the subject lay quietly and 
breathed normally. An HRV analyzer (SSIC, Enjoy Research 
Inc., Taiwan) acquired, stored and processed electrocardio- 
gram signals. Under a sampling rate of 512 Hz, signals were 



recorded using an 8-bit analog-to-digital converter. Stationary 
R-R interval values were resampled and interpolated at a rate 
of 7.11 Hz to produce the continuity in a time domain. Power 
spectral analysis was performed using a nonparametric meth- 
od of fast Fourier transformation (FFT). The direct current 
component was deleted and a Hamming window was used to 
attenuate the leakage effect. The power spectrum was then 
quantified into standard frequency-domain measurements 
defined previously, 34 which consisted of variance (variance of 
RR-interval values), low frequency (LF: 0.04-0.15 Hz), high 
frequency (HF: 0.15-0.40 Hz) and, and the ratio of LF to HF 
(LF/HF). All of the measurements were logarithmically trans- 
formed to correct skewed distribution. 34 Vagal control of HRV 
is represented by HF, whereas both vagal and sympathetic 
control of HRV is jointly represented by LF. The LF/HF ratio 
is considered by some investigators to mirror sympathovagal 
balance or sympathetic modulations. 

Statistical analyses 

SPSS (version 13.0, SPSS, Taipei, Taiwan) statistical soft- 
ware was used for all analyses. Discrete variables in patients 
and controls were compared using chi-square test. Differences 
between continuous variables were evaluated using Student's 
t-test when normally distributed, otherwise the Mann-Whit- 
ney U test was used. Analyses with variances ( ANOVA) with 
post hoc test were used to compare the HRV indices and oth- 
er continuous variables in group comparisons. To evaluate re- 
lations between variables with normal distribution, we used 
Pearson's correlation test. For non-normal distribution Spear- 
man's correlation test was used. The associations between 
HRV measures and age, body mass index (BMI) and habitual 
physical activity were analyzed with product-moment corre- 
lations, whereas point-biserial correlations were used to assess 
relationships with gender. Results of the point-biserial corre- 
lations were identical to those arising from comparisons using 
t tests. Linear regression analyses were used to primarily as- 
sess associations of scores of HAM-A, HAM-D with HRV in- 
dices. To control the confounding effect, we used multiple re- 
gressions on the HRV indices, with HRV-associated factors as 
covariables. All results are two-tailed and a probability value 
p<0.05 was considered statistically significant. 

RESULTS 

Demographics and clinical characteristics 

The groups did not differ significantly on demographic data, 
BMI, systolic/diastolic blood pressure and habitual physical 
activity (Table 1). Three patient groups were comparable re- 
garding age at onset and family history. Group comparisons 
showed significant differences in scores of HAM-A, BAI, 
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Table 1 . Sample characteristics 



Clinical and demographic data 


GAD 


GAD with MD 


MD 


Controls 


Omnibus 
p-value 


Number of participants 


20 


20 


40 


60 




Age, mean±SD, years 


46. 75i 10.95 


47.1 +4.9 


45.88i5.35 


47.45i7.22 


0.75 


Female sex (%) 


16 (80) 


16 (80) 


32 (80) 


48 (80) 


1.0 


Age at onset 


36.95i7.5 


37.1i3.85 


37.1 ±4.6 




0.99 


Family history 


3(15) 


5(25) 


6(15) 


_ 


0.61 


BMI, meaniSD, kg/m 2 


22.59i3.13 


23.71 ±4. 16 


22.96i3.34 


22.41i2.36 


0.4 


Weekly regular exercise, hours 


0.75il.48 


0.2i0.52 


0.98il.l 


0.9il.68 


0.2 


SBP, meaniSD, mm Hg 


122.4±13.8 


121i21.22 


120.3il3.67 


116.28il3.27 


0.3 


DBP, meaniSD, mm Hg 


75.85i 10.64 


73.35ill.04 


74.75i8.84 


72.27i8.22 


0.38 


HAM-A scores, meaniSD 


20.15i6.0 


19.55i5.22 


7.6il.75 


5.42i2.97 


<0.001 


BAI scores, meaniSD 


16.9i9.49 


14.4i8.03 


10.66i2.45 


6.62i5.46 


<0.001 


HAM-D scores, meaniSD 


6.9i2.05 


35.55i8.46 


36.4i6.16 


4.0±2.16 


<0.001 


BDI scores, meaniSD 


8.3il.66 


43 .0+ 10.93 


42.98 i6.48 


4.63i2.96 


<0.001 



Significant comparisons 



GAD vs. controls, GAD with 

MD vs. controls, MD vs. 

controls, GAD vs. MD, GAD 

with MD vs. MD. 
GAD vs. controls, GAD with 

MD vs. controls, MD vs. 

controls, GAD vs. MD, GAD 

with MD vs. MD 
GAD vs. controls, GAD with 

MD vs. controls, MD vs. 

controls, GAD vs. GAD with 

MD, GAD vs. MD. 
GAD vs. controls, GAD with 

MD vs. controls, MD vs. 

controls, GAD vs. GAD with 

MD, GAD vs. MD. 



SD: standard deviation, BMI: body mass index (calculated as weight in kilograms divided by height in meters squared), SBP: systolic blood 
pressure, DBP: diastolic blood pressure, HAM-A: Hamilton Anxiety Rating Scale, BAI: Beck Anxiety Inventory, HAM-D: Hamilton Depres- 
sion Rating Scale, BDI: Beck Depression Inventory, GAD: generalized anxiety disorder, MD: major depression 



HAM-D, and BDI. Post hoc test revealed greater scores of 
HAM-A and BAI for three patient groups vs. controls, GAD 
group vs. MD group and GAD with MD group vs. MD group. 
Greater scores of HAM-D and BDI were shown for three pa- 
tient groups vs. controls, GAD group vs. GAD with MD group 
and GAD group vs. MD group. 

Heart rate variability parameters 

Group comparisons showed significant differences in mean 
R-R intervals (RR), variance, LF and HF (all p values <0.001). 
However, the groups were comparable regarding LF/HF ratio 
(p=0.12). In Figure 1, the GAD with MD group had signifi- 
cantly faster HRs (shorter RR interval) than control group 
(782.451 143.89 milliseconds vs. 923.33±128.66 milliseconds). 
Post hoc test in variance and HF showed lowest values in the 
GAD with MD group [5.94±0.91 ln(ms 2 ) and 3.41±1.05 
ln(ms 2 )], greatest values in the control group [7.29±0.77 
ln(ms 2 ) and 5.48±0.89 ln(ms 2 )], and intermediate values in 



the GAD group [6.65±0.74 ln(ms 2 ) and 4.45±1.34 ln(ms 2 )]. 
Both the GAD group and the GAD with MD group had 
lower LF [4.75±0.9 ln(ms 2 ) and 4.13±1.36 ln(ms 2 )] than the 
control group [5.58±1.13 ln(ms 2 )]. 

Factors associated with HRV 

Associations between HRV measures and those potentially 
confounding variables are summarized in Table 2. Men had 
significantly faster HRs (shorter RR interval) and lower LF 
and LF/HF ratio than women. Participants with higher BMI 
had lower HF. Participants who were habitually more physi- 
cally active had significantly greater LF and HF. The correla- 
tion analysis of individual HRV indices showed that variance 
was positively correlated with mean R-R intervals, LF and HF. 

Association between anxiety/depression severity 
and HRV 

As can be seen in Table 3, subjects with greater scores of 
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HAM-A/HAM-D had significantly faster HRs (shorter RR 
interval). Either HAM-A or HAM-D scores were inversely as- 
sociated with variance, LF and HE However, there was no sig- 
nificant correlation between scores of HAM-A/HAM-D and 
LF/HF ratio. Further adjustment for gender, and BMI did not 



alter the above-mentioned association in a meaningful way 
(Table 4). To analyze the specific influence of anxiety or de- 
pression severity on HRV, we performed a stepwise multiple 
regression with HAM-A/HAM-D scores as predictors of vari- 
ance, LF and HE We found that the HAM-A scores explained 
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Figure 1. Mean R-R Intervals and all measures of HRV for GAD group, GAD with MD group, MD group and controls. Asterisks indicate sig- 
nificant between-groups differences. *p<0.05, **p<0.01, ***p<0.001. Var: total variance [ln(ms 2 )], LF: low frequency power [ln(ms 2 )], HF: 
high frequency power [ln(ms 2 )], LF/HF: ratio of LF to HF [In(ratio)]. 



Table 2. Factors associated with resting HRV indices among all participants 




RR interval 


Var 


LF 


HF 


LF/HF 


Gender (women/men) 1 " 


-0.19* 


-0.19* 


-0.29* 


-0.1 


-0.22* 


Age* 


0.16 


-0.07 


-0.18* 


-0.12 


-0.07 


BMI* 


-0.07 


-0.06 


-0.09 


-0.21* 


0.13 


Physical activity* 


0.12 


0.13 


0.16 


0.14 


0.07 


RR interval 




0.46*** 


0.22* 


0.46*** 


-0.27** 


Var 


0.46*** 




0.83*** 


0.74*** 


0.12 


LF 


0.22* 


0.83*** 




0.62*** 


0.45*** 


HF 


0.46*** 


0.74*** 


0.62*** 




-0.42*** 


LF/HF 


-0.27** 


0.12 


0.45*** 


-0.42*** 





*p<0.05, **p<0.01, ***p<0.001, "fpoint-biserial correlations; first category in parenthesis is the reference group, tproduct-moment correlations. 
BMI: body mass index, Var: total variance [ln(ms 2 )], LF: low frequency power [ln(ms 2 )], HF: high frequency power [ln(ms 2 )], LF/HF: ratio of 
LF to HF [ln(ratio)], HRV: heart rate variability 
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Table 3. The correlations between HRV indices and scores of HAM-Aand HAM-D among all participants, GAD patients, MD patients and 
controls, respectively 

All participants (N=140) GAD (N=40) MD (N=40) Controls (N=60) 

HAM-A HAM-D HAM-A HAM-D HAM-A HAM-D HAM-A HAM-D 



Heart rate variability measures 



RR interval, mean±SD, ms 


-0.27** 


-0.23** 


0.01 


-0.31 


0.14 


0.1 


-0.09 


-0.08 


Var, mean±SD 


-0 37*** 


-0 33*** 


-0.01 


-0.27 


0.13 


0.11 


-0.03 


-0.07 


LF, mean±SD 


-0 3*** 


-0.23** 


-0.18 


-0.25 


0.09 


0.24 


-0.0 


-0.09 


HE meaniSD 


-0.44*** 


-0.38*** 


0.05 


-0.38* 


0.04 


0.16 


-0.09 


-0.21 


LF/HF, meaniSD 


0.15 


0.13 


-0.29 


0.17 


0.1 


0.03 


0.02 


0.07 



*p<0.05, **p<0.01, ***p<0.001. Var: total variance [ln(ms 2 )], LF: low frequency power [ln(ms 2 )], HF: high frequency power [ln(ms 2 )], LF/HF: 
ratio of LF to HF [ln(ratio)], HAM-A: Hamilton Anxiety Rating Scale, HAM-D: Hamilton Depression Rating Scale, GAD: generalized anxi- 
ety disorder, MD: major depression 



Table 4. Multiple regressions of HRV parameters by scores of HAM-A and HAM-D, adjusting for gender, age and BMI factors: unstandard- 
ized regression coefficients 

Unstandardized |3 coefficients (95% confidence interval) 
Unadjusted p Model 1 p Model 2 p Model 3 p 

HAM-A scores 



Mean RR interval 


-0.02 (-0.02-0.01) 


0.001 


-0.02 (• 


■0.02- 


■0.01) 


0.001 


-0.02 (■ 


■0.02-0.01) 


0.001 


-0.02 (-0.03- 


•0.01) 


0.001 


Var 


-3.47 (-4.69-2.25) 


<0.001 


-3.59 (■ 


■4.83- 


■2.35) 


<0.001 


-3.6 (■ 


■4.85-2.35) 


<0.001 


-3.61 (-4.86- 


■2.35) 


<0.001 


LF 


-2.55 (-3.47-1.63) 


<0.001 


-2.78 (■ 


■3.73- 


■1.83) 


<0.001 


-2.85 (• 


■3.81-1.89) 


<0.001 


-2.88 (-3.85- 


•1.9) 


<0.001 


HF 


-2.95 (-3.85-2.06) 


<0.001 


-2.99 (■ 


■3.89- 


■2.08) 


<0.001 


-3.03 (■ 


■3.94-2.11) 


<0.001 


-3.15 (-4.08- 


■2.21) 


<0.001 


LF/HF 


0.44 (-0.73-1.61) 


0.46 


0.46 (• 


■0.74- 


•1.67) 


0.49 


0.46 (• 


•0.75-1.68) 


0.45 


0.45 (-0.78- 


•1.67) 


0.47 


IAM-D scores 


























Mean RR interval 


-0.03 (-0.05-0.01) 


0.009 


-0.03 (■ 


■0.05- 


■0.01) 


0.007 


-0.03 (• 


■0.05-0.01) 


0.009 


-0.03 (-0.05- 


■0.01) 


0.016 


Var 


-6.2 (-9.12-3.28) 


<0.001 


-6.5 (• 


■9.47- 


■3.52) 


<0.001 


-6.56 (■ 


•9.54-3.58) 


<0.001 


-6.3 (-9.26- 


■3.34) 


<0.001 


LF 


-3.34 (-5.59-1.08) 


0.004 


-3.73 (■ 


■6.08- 


■1.37) 


0.002 


-3.94 (■ 


■6.32-1.56) 


0.001 


-3.64 (-6.02- 


■1.26) 


0.003 


HF 


-4.72 (-6.94-2.49) 


<0.001 


-4.8 (■ 


■7.04- 


■2.56) 


<0.001 


-4.95 (• 


■7.2-2.7) 


<0.001 


-4.59 (-6.89- 


■2.28) 


<0.001 


LF/HF 


1.68 (-1.01-4.37) 


0.22 


1.69 (■ 


■1.08- 


■4.45) 


0.23 


1.65 (• 


■1.12-4.42) 


0.24 


1.34 (-1.42- 


■4.1) 


0.34 



Model 1: controlling for gender, Model 2: controlling for model 1+age, Model 3: controlling for model 1+age+BMI, Var: total variance 
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Anxiety Rating Scale, HAM-D: Hamilton Depression Rating Scale, BMI: body mass index 



significantly the variation of HF [r 2 =0.234, F(l, 138)=42.22, 
p<0.001, b=-2.95] (Figure 2A), variance [r^O.187, F(l, 138)= 
31.64, p<0.001, b=-3.47] and LF [r^O.18, F(l, 138)=30.27, 
p<0.001, b=-2.55]. To a lesser extent, The HAM-D scores ex- 
plained significantly the variation of HF [r 2 =0.113, F(l, 
138)=17.6, p<0.001, b=-4.72] (Figure 2B), variance [r^O.113, 
F(l, 138)=17.58, p<0.001, b=-6.2] andLF [r^O.06, F(l, 138)= 
8.54, p=0.004, b=-3.34]. When separately analyzing the corre- 
lations between clinical symptoms and HRV indices among 
GAD patients, MD patients and controls, we only found 
HAM-D scores negatively correlated with HF among GAD 
patients. However, further adjustment for gender, age, BMI 
and physical activity largely reduced the (3 coefficients and 
rendered the correlation non-significant (p=0.075). 



DISCUSSION 

The main results of our study are summed as follows. First, 
we found physically healthy, unmedicated GAD patients had 
lower HRV than healthy controls. This result is consistent 
with previous investigations 17,20 suggesting that reduction in 
HRV is a psychophysiological marker of GAD. However, the 
largest study published so far reported that GAD patients did 
not differ from controls in HRV and that lower HRV in cur- 
rent GAD was derived from the effect of antidepressants rath- 
er than the diagnosis itself. 22 We believe that our result is reli- 
able based on the following reasons. 1) Although the sample 
size of the study undertaken by Licht and colleagues is impres- 
sive, a number of methodological issues may account for their 
findings. For example, the authors only used time-domain 
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measures of HRV. It is possible that very subtle autonomic 
changes accompany GAD and a time-domain analysis of HRV 
is insufficient to detect such small differences. Frequency-do- 
main analysis of HRV is a sophisticated tool for the detection 
of the cardiac autonomic regulation and can add information 
on quantification of parasympathetic and sympathetic ner- 
vous system function over time-domain analysis of HRV 1,35 
2) We have excluded subjects with psychiatric and physical 
co -morbidities that could potentially confound the associa- 
tion between GAD and cardiac autonomic functions. Partici- 
pants in the GAD group were interviewed with the modified 
Chinese Version of SADSL 29 to rule out any psychiatric co- 
morbidity. Thus a false -positive result due to inclusion of de- 
pressive disorders, other anxiety disorders or substance use 
disorder in the GAD group is presently unlikely. 3) We have 
also controlled other confounding factors that may suppress 
or magnify the true effects of GAD on HRV, including medi- 
cation, age, smoking, BMI and physical activity levels. 15 For 
example, all our participants do not smoke or have a history 
of smoking. Current smoking clearly depresses HRV and, 
even among those who have recently quit, HRV remains low- 
er compared with that of normal nonsmokers. 36 Therefore, 
smoking must be taken into account in any study of the effect 
of GAD on HRV because a substantial proportion of GAD 
patients smoke or have a history of smoking. 37 4) Ethnic strat- 
ification among study samples may lead to resetting popula- 
tion HRV patterns. 38,39 For example, evidence indicated that 
Asian ethnicity served as significant predictors of healthy 
HRV as compared to European-Americans in a non-medical 
sample. 39 Studying a genetically heterogeneous population 
might produce a false-positive or false-negative result by 
chance rather than reveal a direct relation. However, all our 
subjects were unrelated Han Chinese subjects drawn from a 
population pool in Taiwan that is known to be genetically ho- 
mogeneous. 41 All of the biological grandparents of our re- 
cruited subjects were of Han Chinese ancestry. Therefore, it is 
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Figure 2. A stepwise multiple regression 
showed that the HAM-A scores explains 
23.4% of HF variation (A), and the HAM- 
D scores explains 11 .3% of HF variation 
(B). HF: high frequency power [ln(ms 2 )], 
HAM-A: Hamilton Anxiety Rating Scale, 
HAM-D: Hamilton Depression Rating 
Scale. 



less likely that ethnic stratification bias produced a false-posi- 
tive result in our study. 

Second, correlation analysis showed a strongly positive as- 
sociation between variance and HF-HRV (r=0.73, p<0.001), 
suggesting that an overall reduction in HRV in GAD patients 
was at least partly derived from suppression of parasympa- 
thetic input as indicated by low HF-HRV. The following two 
findings further complement our result regarding lower HF- 
HRV in GAD. 1) Patients with more severe anxiety symptoms 
tend to have lower HF-HRV than those with less severe anxi- 
ety symptoms. This finding is consistent with prior work re- 
porting that GAD and its cardinal feature (worry), are associ- 
ated with lower cardiac vagal control. 17 2) Given the large 
contribution of the PNS to resting HR, findings for resting 
HF-HRV in any single study will theoretically parallel the 
findings for resting HR. This theoretical perspective seems to 
be supported by the fact that resting HR in GAD patients ele- 
vated on a trend level (Figure 1). However, it would be pre- 
mature to make a conclusion regarding resting HR in GAD 
based up the preliminary result reported herein. In order to 
conclude GAD-related alterations in resting HR, larger sam- 
ples and/or meta-analyses are needed. Taken together, our 
findings and can be interpreted as evidence that the basal au- 
tonomic state of GAD is characterized by decreased parasym- 
pathetic tone. The following theoretical perspectives support 
our findings. 1) Polyvagal theory proposed by Porges high- 
lights the importance of the vagal pathway in attention, emo- 
tion expression, social bonding and flexible adjustment to en- 
vironmental demands, 41 and all of them are compromised in 
GAD patients. 42 Vagal influences to the heart serve to dampen 
the sympathetic reactions to stress and to promote calm be- 
havioral states and self-regulation. 14 Without this protective 
function of vagal tone, subjects may become vulnerable to 
anxious apprehension and worry, involving pre-attentive bi- 
ases to threat information, and rigid and inflexible response 
patterns, which characterized GAD. 42 2) Another explanation 
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proposed for the link between GAD and decreased parasym- 
pathetic tone is the inability to disengage threat detection, 
which serves to perpetuate hyperarousal and worry, even 
when no real threat exists, 43 which may cause a chronic with- 
drawal of PNS activity and long-term reductions in HRV, 
which subsequently increasing the risk for CVD. In consis- 
tent with this explanation, two prospective follow-up studies 
on patients with stable CAD reported GAD as having a dis- 
tinctly negative impact on subsequent cardiovascular 
events. 6,7 

Third, significant reduced LF-HRV in our GAD patients 
did not cohere with what we would have expected from other 
cardiac autonomic results, i.e. elevated LF-HRV 17 In addition, 
we failed to find a LF/HF ratio difference between GAD pa- 
tients and controls, although it is hypothesized that reduced 
vagal modulation is accompanied by a subsequent displace- 
ment of the sympathovagal balance (as indexed by the LF/HF 
ratio) in favor of sympathetic modulation. These unexpected 
findings warrant comment. 1) The traditional interpretations 
of the HRV measures used in our study are that HF power es- 
timates vagal tone, while LF power reflects both vagal and 
sympathetic influences. However, it has also been reported 
that when LF power is assessed in the supine position, ad- 
ministration of atropine (a potent inhibitor of parasympa- 
thetic muscarinic receptors) eliminates most of the LF region 
of the power spectrum. 44 This does not occur when LF power 
is assessed in the sitting position, and suggests that resting LF 
power in our study may primarily reflect vagal influences. 45 
Consistent with this speculation is evidence that there is a 0.62 
correlation between LF-HRV and HF-HRV (Table 2). 2) The 
LF/HF ratio was calculated from the absolute values of the LF 
and HF power for each subject and was considered to be an 
indicator of the sympathovagal balance or sympathetic mod- 
ulations. However, there is mounting evidence against LF/HF 
ratio representing sympathovagal balance. For example, Goed- 
hart et al. 46 reported that the LF/HF ratio did not show the 
expected correlation to the preejection period, an established 
measure of cardiac sympathetic control. More specifically, 
Goldstein et al. 47 considered LF power to be an index not of 
cardiac sympathetic tone but of baroreflex function during 
supine rest. Additional studies are recommended (e.g., cardiac 
noradrenaline spillover) for adequately assessing sympathetic 
nerve activity 

Finally, we found that HF-HRV was more reduced in GAD 
with MD group as compared to GAD group, and that these 
reductions were not due to anxiety severity - both groups rat- 
ed similarly on levels of anxiety severity - or other potential 
confounding variables such as age, gender, BMI and physical 
activity. This finding is reinforced by the inverse correlation 
between HF-HRV and depression severity (Figure 2B), which 



coheres with what we would have expected from our prior 
study result. 26 Our data provides an important extension to 
prior work on GAD, 17,20,21 indicating that comorbid depres- 
sion is associated with the greater reductions in CVC. In con- 
trast to our findings, Hofmann et al. 48 observed that GAD pa- 
tients with comorbid MD had greater HF-HRV values than 
did those without MD. However, the authors use medicated 
patient sample, and therefore the influence of medication 
cannot be ruled out. Moreover, there was no healthy control 
group in that study. Lack of such a comparison with controls 
would have weakened the interpretability of their findings. 
We proposed two explanations for why HF-HRV may be re- 
duced in GAD and in those with comorbid MD in particular. 
One potential explanation is the cumulative effect, i.e. com- 
bining the psychophysiological impact of the two disorders. 
Previous research has shown that both GA D 17 ' 21 and MD 15,26 
are associated with decreased CVC. A plausible scenario may 
be that patients with comorbid GAD and MD have greater 
reduction in CVC than those with only one diagnosis. As can 
be seen in Figure 1, the GAD with MD group had the lowest 
values of HF-HRV among three patient groups indeed. Also, 
the cumulative effect may be partly responsible for the signifi- 
cantly elevated HRs in the GAD with MD group (Figure 1). 
Another explanation is that the comorbidity between the two 
disorders may recognize a unique group of patients. 49 Indeed, 
the nosological rules of DMS-IV reflect the uncertainty as to 
whether GAD and MD represent two separate diagnostic en- 
tities. For example, the symptoms of GAD do not occur ex- 
clusively during a major depressive episode while one can 
only made a diagnosis of GAD in individuals with MD if the 
GAD symptoms also occur outside a major depressive epi- 
sode. Our data provide further evidence for the possibility 
that individuals with comorbid GAD and MD can be distin- 
guished based on psychophysiological correlates (for example, 
HF-HRV) from individuals with GAD but without a second- 
ary diagnosis of MD. One might argue that the differences in 
HRV indices between the two groups reflected difference in 
depression severity because HRV indices were significantly 
correlated with depression severity in all of our participants 
(Table 3). However, separate analysis of 40 GAD patients and 
further adjustment for HRV-associated covariables showed 
that HRV indices were not associated with depression severi- 
ty Taken as a whole, the between-group HRV differences did 
not proxy differences in severity but reflected psychophysio- 
logical differences between non-comorbid GAD and GAD 
comorbid with depression. Moreover, it is worth mentioning 
that Frasure- Smith et al. implicated GAD alone, MD alone 
and concurrent GAD and MD as all increasing the possibility 
of a poor cardiac outcome among CAD patients. However, 
patients with comorbid GAD and MDD were not at greater 
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risk for major adverse cardiac events than those with only 1 
factor. 7 Therefore, our result must be considered preliminary. 
Further investigation of the link between GAD with comor- 
bid depression, HRV and CVD is warranted. 

As suggested above, our findings highlight low cardiac va- 
gal tone in GAD patients and in those with comorbid MD in 
particular and should serve to remind clinicians to pay atten- 
tion to their increased risk of CVD. For example, an auto- 
nomic function examination such as HRV analysis can be 
done to provide a rapid screening of systemic autonomic dis- 
turbance. GAD patients with comorbid depression may ben- 
efit better from cardiovascular risk reduction strategies such 
as exercise, smoking cessation, or undergoing treatment to 
restore the autonomic function. For instance, mindfulness 
meditation has been reported to significantly increase cardiac 
parasympathetic activity 50 and reduce symptoms of GAD. 51 
Moreover, cognitive behavior therapy-based psychotherapy 
applied in a forest environment was helpful in the achieve- 
ment of depression remission while increased parasympa- 
thetic nerve tone. 52 Larger randomized control trials in GAD 
patients with comorbid depression to validate these study re- 
sults are warranted. 

Several limitations should be considered in the present 
study. 1) Interpretation of LF-HRV as sympathetic activity 
should be made with caution, as the level of modulation of 
LF-HRV by the sympathetic branch is still debated. 2) Our 
study failed to provide the information concerning the men- 
strual cycle of our female participants; however, HRV was 
shown to fluctuate during different phases of female men- 
strual cycle. 53 3) Our study did not analyze time-domain 
measures of HRV, which might provide more clues regarding 
the different findings between our study and the study by Li- 
cht and colleagues. 22 4) A recent meta-analysis emphasizes 
the importance of nonlinear techniques to study cardiac auto- 
nomic function and appear to be superior or of additional val- 
ue to the traditional time and frequency domain measures. 54 
Thus, further replication with nonlinear measures is clearly 
desirable. 

Acknowledgments 

This study was supported by Tri-Service General Hospital Grant TSGH- 
C98-91 and TSGH-C102-123 (H.A.C). The authors thank psychologist 
Tsuey-Yen Yeh for her assistance in the interviewing and Miss Zhang- Yu 
Wang for her assistance in preparing this manuscript. 

REFERENCES 

1. Cohen H, Benjamin J. Power spectrum analysis and cardiovascular 
morbidity in anxiety disorders. Auton Neurosci 2006;128:1-8. 

2. Liao D, Cai J, Rosamond WD, Barnes RW, Hutchinson RG, Whitsel 
EA, et al. Cardiac autonomic function and incident coronary heart dis- 
ease: a population-based case-cohort study. The ARIC Study. Athero- 
sclerosis Risk in Communities Study. Am J Epidemiol 1997;145:696- 
706. 



3. Tsuji H, Larson MG, Venditti FJ Jr, Manders ES, Evans JC, Feldman CL, 
et al. Impact of reduced heart rate variability on risk for cardiac events. 
The Framingham Heart Study. Circulation 1996;94:2850-2855. 

4. Coryell W, Noyes R, Clancy J. Excess mortality in panic disorder. A 
comparison with primary unipolar depression. Arch Gen Psychiatry 
1982;39:701-703. 

5. Coryell W, Noyes R Jr, House JD. Mortality among outpatients with 
anxiety disorders. Am J Psychiatry 1986;143:508-510. 

6. Martens ET, de Jonge P, Na B, Cohen BE, Lett H, Whooley MA. Scared 
to death? Generalized anxiety disorder and cardiovascular events in pa- 
tients with stable coronary heart disease: The Heart and Soul Study. 
Arch Gen Psychiatry 2010;67:750-758. 

7. Frasure-Smith N, Lesperance F. Depression and anxiety as predictors of 
2-year cardiac events in patients with stable coronary artery disease. 
Arch Gen Psychiatry 2008;65:62-71. 

8. Elghozi JL, Girard A, Laude D. Effects of drugs on the autonomic con- 
trol of short-term heart rate variability. Auton Neurosci 2001;90:116- 
121. 

9. Hoehn-Saric R, McLeod DR, Zimmerli WD. Somatic manifestations in 
women with generalized anxiety disorder. Psychophysiological re- 
sponses to psychological stress. Arch Gen Psychiatry 1989;46:1113- 
1119. 

10. Hofmann SG, Moscovitch DA, Litz BT, Kim HJ, Davis LL, Pizzagalli 
DA. The worried mind: autonomic and prefrontal activation during 
worrying. Emotion 2005;5:464-475. 

11. Roth WT, Doberenz S, Dietel A, Conrad A, Mueller A, Wollburg E, et 
al. Sympathetic activation in broadly defined generalized anxiety disor- 
der. J Psychiatr Res 2008;42:205-212. 

12. Fisher AJ, Granger DA, Newman MG. Sympathetic arousal moderates 
self-reported physiological arousal symptoms at baseline and physio- 
logical flexibility in response to a stressor in generalized anxiety disor- 
der. Biol Psychol 2010;83:191-200. 

13. Kalk NJ, Nutt DJ, Lingford-Hughes AR. The role of central noradrener- 
gic dysregulation in anxiety disorders: evidence from clinical studies. J 
Psychopharmacol 2011;25:3-16. 

14. Porges SW. The polyvagal theory: phylogenetic substrates of a social 
nervous system. Int J Psychophysiol 2001;42:123-146. 

15. Rottenberg J. Cardiac vagal control in depression: a critical analysis. 
Biol Psychol 2007;74:200-211. 

16. Friedman BH. An autonomic flexibility-neurovisceral integration 
model of anxiety and cardiac vagal tone. Biol Psychol 2007;74:185-199. 

17. Thayer JF, Friedman BH, Borkovec TD. Autonomic characteristics of 
generalized anxiety disorder and worry. Biol Psychiatry 1996;39:255- 
266. 

18. Heart rate variability: standards of measurement, physiological inter- 
pretation and clinical use. Task Force of the European Society of Cardi- 
ology and the North American Society of Pacing and Electrophysiolo- 
gy. Circulation 1996;93:1043-1065. 

19. van Zyl LT, Hasegawa T, Nagata K. Effects of antidepressant treatment 
on heart rate variability in major depression: a quantitative review. Bio- 
psychosoc Med 2008;2:12. 

20. Kollai M, Kollai B. Cardiac vagal tone in generalised anxiety disorder. 
Br J Psychiatry 1992;161:831-835. 

21. Pittig A, Arch J J, Lam CW, Craske MG. Heart rate and heart rate vari- 
ability in panic, social anxiety, obsessive-compulsive, and generalized 
anxiety disorders at baseline and in response to relaxation and hyper- 
ventilation. Int J Psychophysiol 2013;87:19-27. 

22. Licht CM, de Geus EJ, van Dyck R, Penninx BW. Association between 
anxiety disorders and heart rate variability in The Netherlands Study of 
Depression and Anxiety (NESDA). Psychosom Med 2009;71:508-518. 

23. Breslau N, Schultz L, Peterson E. Sex differences in depression: a role 
for preexisting anxiety. Psychiatry Res 1995;58:1-12. 

24. Kessler RC, Nelson CB, McGonagle KA, Liu J, Swartz M, Blazer DG. 
Comorbidity of DSM-III-R major depressive disorder in the general 
population: results from the US National Comorbidity Survey. Br J Psy- 



334 Psychiatry Investig 2013;10:326-335 



HA Chang et al. 



chiatry Suppl 1996;30:17-30. 

25. Kessler RC, DuPont RL, Berglund P, Wittchen HU. Impairment in pure 
and comorbid generalized anxiety disorder and major depression at 12 
months in two national surveys. Am J Psychiatry 1999;156:1915-1923. 

26. Chang HA, Chang CC, Chen CL, Kuo TBJ, Lu RB, Huang SY. Major 
depression is associated with cardiac autonomic dysregulation. Acta 
Neuropsychiatrica 2012;24:318-327. 

27. Song BA, Yoo SY, Kang HY, Byeon SH, Shin SH, Hwang EJ, et al. Post- 
traumatic stress disorder, depression, and heart-rate variability among 
north Korean defectors. Psychiatry Investig 201 1;8:297-304. 

28. Chang HA, Chang CC, Tzeng NS, Kuo TB, Lu RB, Huang SY. De- 
creased cardiac vagal control in drug-naive patients with posttraumatic 
stress disorder. Psychiatry Investig 2013;10:121-130. 

29. Endicott J, Spitzer RL. A diagnostic interview: the schedule for affective 
disorders and schizophrenia. Arch Gen Psychiatry 1978;35:837-844. 

30. Huang SY, Lin WW, Ko HC, Lee JF, Wang TJ, Chou YH, et al. Possible 
interaction of alcohol dehydrogenase and aldehyde dehydrogenase 
genes with the dopamine D2 receptor gene in anxiety-depressive alco- 
hol dependence. Alcohol Clin Exp Res 2004;28:374-384. 

31. Beck AT, Steer RA. Beck Anxiety Inventory Manual. San Antonio: The 
Psychological Corporation Harcourt; 1990. 

32. Hamilton M. The assessment of anxiety states by rating. Br J Med Psy- 
chol 1959;32:50-55. 

33. Beck AT, Ward CH, Mendelson M, Mock J, Erbaugh J. An inventory 
for measuring depression. Arch Gen Psychiatry 1961;4:561-571. 

34. Kuo TB, Lin T, Yang CC, Li CL, Chen CF, Chou P. Effect of aging on 
gender differences in neural control of heart rate. Am J Physiol 1999; 
277:H2233-H2239. 

35. Akselrod S, Gordon D, Ubel FA, Shannon DC, Berger AC, Cohen RJ. 
Power spectrum analysis of heart rate fluctuation: a quantitative probe 
of beat-to-beat cardiovascular control. Science 1981;213:220-222. 

36. Stein PK, Rottman JN, Kleiger RE. Effect of 21 mg transdermal nico- 
tine patches and smoking cessation on heart rate variability. Am J Car- 
diol 1996;77:701-705. 

37. Moylan S, Jacka FN, Pasco JA, Berk M. Cigarette smoking, nicotine de- 
pendence and anxiety disorders: a systematic review of population- 
based, epidemiological studies. BMC Med 2012;10:123. 

38. Li Z, Snieder H, Su S, Ding X, Thayer JF, Treiber FA, et al. A longitudi- 
nal study in youth of heart rate variability at rest and in response to 
stress. Int J Psychophysiol 2009;73:212-217. 

39. Martin LA, Doster JA, Critelli JW, Lambert PL, Purdum M, Powers C, 
et al. Ethnicity and Type D personality as predictors of heart rate vari- 
ability. Int J Psychophysiol 2010;76:1 18-121. 

40. Chang CC, Lu RB, Chen CL, Chu CM, Chang HA, Huang CC, et al. 
Lack of association between the norepinephrine transporter gene and 
major depression in a Han Chinese population. J Psychiatry Neurosci 
2007;32:121-128. 



41. Porges SW. Orienting in a defensive world: mammalian modifications 
of our evolutionary heritage. A Polyvagal Theory. Psychophysiology 
1995;32:301-318. 

42. Bradley BP, Mogg K, White J, Groom C, de Bono J. Attentional bias for 
emotional faces in generalized anxiety disorder. Br J Clin Psychol 1999; 
38:267-278. 

43. Thayer JF, Lane RD. A model of neurovisceral integration in emotion 
regulation and dysregulation. J Affect Disord 2000;61:201-216. 

44. Pomeranz B, Macaulay RJ, Caudill MA, Kutz I, Adam D, Gordon D, et 
al. Assessment of autonomic function in humans by heart rate spectral 
analysis. Am J Physiol 1985;248:H151-H153. 

45. Sloan RP, McCreath H, Tracey KJ, Sidney S, Liu K, Seeman T. RR inter- 
val variability is inversely related to inflammatory markers: the CAR- 
DIA study. Mol Med 2007;13:178-184. 

46. Goedhart AD, Willemsen G, Houtveen JH, Boomsma DI, De Geus EJ. 
Comparing low frequency heart rate variability and preejection period: 
two sides of a different coin. Psychophysiology 2008;45:1086-1090. 

47. Goldstein DS, Bentho O, Park MY, Sharabi Y. Low-frequency power of 
heart rate variability is not a measure of cardiac sympathetic tone but 
may be a measure of modulation of cardiac autonomic outflows by 
baroreflexes. Exp Physiol 2011;96:1255-1261. 

48. Hofmann SG, Schulz SM, Heering S, Muench F, Bufka LF. Psychophys- 
iological correlates of generalized anxiety disorder with or without co- 
morbid depression. Int J Psychophysiol 2010;78:35-41. 

49. Mineka S, Watson D, Clark LA. Comorbidity of anxiety and unipolar 
mood disorders. Annu Rev Psychol 1998;49:377-412. 

50. Takahashi X Murata T, Hamada T, Omori M, Kosaka H, Kikuchi M, et 
al. Changes in EEG and autonomic nervous activity during meditation 
and their association with personality traits. Int J Psychophysiol 2005; 
55:199-207. 

5 1 . Wong SY, Mak WW, Cheung EY, Ling CY, Lui WW, Tang WK, et al. A 
randomized, controlled clinical trial: the effect of mindfuhiess-based 
cognitive therapy on generalized anxiety disorder among Chinese com- 
munity patients: protocol for a randomized trial. BMC Psychiatry 201 1; 
11:187. 

52. Kim W Lim SK, Chung EJ, Woo JM. The effect of cognitive behavior 
therapy-based psychotherapy applied in a forest environment on physi- 
ological changes and remission of major depressive disorder. Psychiatry 
Investig 2009;6:245-254. 

53. Sato N, Miyake S, Akatsu J, Kumashiro M. Power spectral analysis of 
heart rate variability in healthy young women during the normal men- 
strual cycle. Psychosom Med 1995;57:331-335. 

54. Kemp AH, Quintana DS, Gray MA, Felmingham KL, Brown K, Gatt 
JM. Impact of depression and antidepressant treatment on heart rate 
variability: a review and meta-analysis. Biol Psychiatry 2010;67:1067- 
1074. 



www.psychiatryinvestigation.org 335 



